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I. INTRODUCTION 



Despite the remarkable success of the gauge sector of the Standard Model (SM), there still exist 
some problems in the Higgs and Yukawa sectors. The Yukawa matrices are responsible for the 
masses and mixings of matter fermions: quarks and leptons. The Yukawa sector in the SM can 
give experimentally consistent masses and mixings, because it contains more free parameters than 
the number of observables in general. There is no predictivity in the Yukawa sector because of this 
redundancy of the parameters. One of the ideas to overcome this issue is to introduce a family 
symmetry (flavor symmetry), which is the symmetry between generations. In this paper we consider 
a concrete model which is symmetric under the binary dihedral group Qq [l|, 0]. 

On the other hand, in the Higgs sector, the most important problem is that the Higgs boson 
has not been experimentally discovered yet. Discovery of the Higgs boson is expected at the Large 
Hadron Collider (LHC). In the SM, the Higgs mass is quadratic divergent. This problem is solved 
by introducing Supersymmetry (SUSY) at 0(1) TeV. The Minimal Supersymmetric Standard Model 
(MSSM) has low energy SUSY. In general it contains gauge symmetric, lepton and baryon number 
violating operators 



Wu 



-XijkLiLjEl + \' ijh .L,( ):!){. + -X'l^U-D^Dl + ^H u Li 



in addition to the usual Yukawa couplings and /z-term. The asymmetric properties Ay& 



(1) 

-X jik and 

Kjk = ~Kkj mean t na t 9 + 27 + 9 + 3 = 48 (complex) parameters are included in this interactions. 
These couplings generate unacceptable processes such as Lepton Flavor Violating (LFV) processes 
and proton decay. The conservation of R-parity 



R 



-1) 



3B+L+2s 



(2) 



where B, L and s denote baryon, lepton number and spin of the particles, respectively, is one 
possibility to forbid the couplings Eq.([T]). From the definition, R-parity is +1 for all SM particles 
and —1 for all their superpartners. However, R-parity is not the only possible choice to forbid 
the interactions. Matter- or lepton- and baryon-parity [13] can be also a possibility. On the other 
hand, without R-parity, these coupling constants have to be strictly constrained not to conflict with 



experimental data. Constraints on the R-parity violating couplings 
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rave been obtained by many 
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proton decay [16|, l33_|, [34 1 , and so on 



Family symmetries also constrain the form of R-parity interactions 
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23 L neutral meson 



as well as the 




-parity 

Yukawa matrices. In the model that we consider [l|, |2j, the Qq family symmetry reduces the 45 
trilinear couplings to three: A, X[ and X' 2 . The baryon number violating couplings X'-j k UfDjDl are 
forbidden by the symmetry in our model, so it is guaranteed by the symmetry that the R-parity 
violating operators do not induce proton decay. 
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In this paper, we study the phenomenology of the three R-parity violating interactions in the 
model with Q e family symmetry l|, Q|. First, we obtain upper bounds on three coupling constants 
A and X[ 2 from the experimental constraints. Next we focus on LFV processes induced by A. The 
XLLE C coupling generates the LFV decays £~ — > £~£~£^(m,i,j,k denote the flavor of the charged 
lepton) at tree level, and the Branching Ratios (BR) of the decay processes are proportional to A 4 . 
Therefore, the ratios of these processes are independent of A and can be predicted unambiguously to 
be BR(t — > eee)/ BR{t — > nnn) ~ 4m^/mJ. It reflects the properties of the family symmetry. We 
introduce the Q$ symmetric model in the next section, and derive the predictions in the sect. III. 

II. THE MODEL 

A. Group Theory of Qq 

The binary dihedral group Qn(N = 2,4,6, ...) is a finite subgroup of SU(2) and defined by the 
following set of 2N elements 

Q N = {\ t A i A 2 t ...,A N - l ,B,AB,... i A N - 1 B} , (3) 

where two dimensional representation of matrix A and B is given by 

A=( = B={ { V (4) 

\ —sm(p N cos0Ar J JM \ — i J 

Since we consider a supersymmetric model with Q§ family symmetry, we show the multiplication 
rules only for the case of N = 6. 

Qe group contains 2 two-dimensional irreducible representations (irreps), 2\,2 2 and 4 one- 
dimensional ones l+,o, l+,2> l-,3) where 2i is pseudo real and 2 2 is real representation. In 
the notation of l± n (n = 0,1,2,3), ± stands for the change of sign under the transformation by 
matrix A, and n the factor exp(m7r/2) by B. So l +) o and l +j 2 are real representations, while l_ i 
and 1_ 3 are complex conjugate to each other. Their group multiplication rules are given as follows 
UBS': 

1 

(5) 



l+,o + 1 .2 2 2 

-x x y 2 - x 2 yi \ (6) 
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,1 
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\ V2 



(xm - x 2 yi) (xm + x 2 y 2 ) 



xm ~ x 2 y 2 
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Oi 1 
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a 2 J 
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{aib x 



+ 1+.2 + 
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(xia 2 + x 2 ai) (xiai - x 2 a 2 ) 
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Xiai + rc 2 a 2 
Xia 2 — x 2 ai 



(8) 



In what follows we construct a concrete model with Q e family symmetry by using of the multiplication 
rules above. 

B. Qq assignment and superpotential 

We show the Qq assignment of the quark, lepton and Higgs chiral supermultiplets below, where 
Qi, Q3, Li, L3 and Hf, H%, Hf, H% stand for SU (2)l doublets supermultiplets for quarks, leptons and 
Higgs bosons, respectively. Similarly, SU(2)l singlet supermultiplets for quarks, charged leptons and 
neutrinos are denoted by Uf, U£, Dj, D 3 , Ej, E% and Nf, N£. The generation indices J, J, ... = (1, 2) 
are applied to the Q 6 doublet, and ... = (1, 2, 3) to three generations throughout the paper. Y is 
gauge singlet Higgs supermultiplet to give neutrino mass by seesaw mechanism. 

We give the Qq assignment to each field as follows: 



2i 


Qi, 








2 2 


uj, 


Dj, 


U, 


Ej, Nj, 


l+,o 




El 








Qs, 


Y, 






1 -.1 


Us, 


D c 3 , 




Tjd 


l-,3 


N c 3 . 









(9) 



In a model without R-parity conservation, there is no distinction between lepton doublet and down 
type Higgs doublet, because both have the same gauge quantum numbers. We have written these 
fields as Li, Hf, and physical lepton doublet and down type Higgs doublet will be written as linear 
combination of these. At first we write down the superpotential in the fields with hat (Lj, Hf), and 
after that we rewrite it in the physical fields without hat (Lj, Hf). 

Under the field assignment above, we can write down the most general, renormalizable Q§ invariant 
superpotential W (without R-parity conservation): 



W = W Q + W L + W^ + Wp, 



(10) 
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where 



W Q = Wu + W D , 

W v = Y a u Q 3 U^ + Y b u Qj {a^jjinBTj -YfQ 3 Uf (ta 2 ) u H u j + Y C U Q, (a 1 ) u U c jHl 
W D = Y a d Q 3 D^ + Y b d Q I (a 1 ) IJ D^-Y b lQ 3 D c I ( l a 2 ) IJ Hj^ 

(11) 

for the quark sector 1 , and 
W L = W E + W N , 

W E = Y b e LjEIHf + Y b 1L 3 EjHf + Y c e fiJKL I E^ J H d K , 

W N = Y:L 3 N C 3 H% + Y b "L 3 NjHj + YZf JJK LiN5H% + ^M N N<N< + \ N YN*N C 3 , 

—fin = /221 = /122 = /212 = 1 (12) 

for the lepton sector. The R-parity violating trilinear coupling Wp which respects the family sym- 
metry is given by 

W n = XL.LjE'j + \' 1 L I (ia 2 ) IJ Q 3 D<j + ^L^uQjD^. (13) 

The family symmetry constrains not only the form of the Yukawa sector but also that of the R-parity 
violating couplings. Proton decay caused by dimension-four operators X'-^U^DjD^ are prevented 
since the baryon number violating interactions are forbidden by the symmetry jj]. Therefore R- 
parity violating trilinear terms Wn contain only three couplings A and A' 12 . It should be compared 
with the MSSM case in which it contains 45 (complex) trilinear couplings. We will analyze the 
phenomenology of the couplings Eq. ( jTBl) in the section III. 

In the following analysis, we assume that any couplings appearing in the superpotential Eq. ffTOl) 
are real, and CP symmetry is violated spontaneously by vacuum expectation values (VEVs) of the 
Higgs bosons. The //-term W^ which includes both R-parity conserving and violating bilinear terms 
is discussed in the next subsection. 

C. Bilinear terms 

In this subsection, we discuss bilinear terms and define the physical fields Lj and Hj. Since the 
lepton and down type Higgs doublet have the same gauge quantum numbers, it is natural that both 



The superpotential for the up quark sector Wu does not have hat ("), because it contains neither Li nor Hf. 
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Li and Hf get VEVs. We assume that the fields get complex VEVs 

W> = (#2> = (# 3 ") = Jf* e **> ( H £) = 7f V 3^, (14) 

{Hi) = {Hf) = \v d e ie \ (U) = (U) = \v L e id \ (15) 

where the VEVs of Hf and Lj have the same phase. Also, the Qq singlet lepton doublet L 3 does 
not have non-zero VEV by definition. To ensure this VEV structure, the //-term should have the 
symmetry 

HI <-> HI, Hf <-> Hi U <-> £2, (16) 

and it is written as 2 

# M = fifH^Hf + riBZH* + AziHl + H u 2 )H d , + ^(Hf + ^2) + ^H^ijH* 

+ ^iZyLj + /4i#s (^i + La) + /^tf/V)^./. (17) 

We define the physical Higgs fields as the fields whose VEVs break SU(2)l X U(l)y symmetry. So 
the physical Higgs fields Hf and lepton doublets Lj are defined by the linear combination 

H i = -fH, + -jLr, L l = -jH 1 - -^Lj, (18) 
u D u D u D u D 

where v d D = vA^J-P v\. Therefore VEVs of Hf are 

(Hf) = V fe^ d , (19) 

which break the electroweak symmetry, while Lj do not get VEVs. 

The /i-term Eq. ffTT|) should be rewritten in the new Higgs fields Hf and lepton doublet Lj as 
follows: 

W, = m coa&HZHf + + ^ 3 (#* + H%)Hl + // 31 cos£ 3 i# 3 u (#i + tf 2 d ) + /i 12 cos £ 12 H? (a 1 ) uH* 

+ m sm^HfLj + M 3i sin^i^ 3 U (Li + £2) + a* 12 sin ^ 2 Hf(a l ) u Lj, (20) 

where the mixing angles of the Higgs fields and lepton doublets are defined as 

//isin^i = -3 , etc. (21) 



2 This form of the /z-term is given by introducing additional gauge singlet Higgs fields and extra discrete symmetry 
{HQ. In the present paper, we just assume the form of Eq. (fTT)) . 
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with /ix = A/(/if) 2 + (/if) 2 etc. 

These mi xing ter ms of neutralinos and neutrinos generate neutrino masses proportional to sin 2 £ 
at tree level 



21 



22]. 



There are particularly strong constraints on the mixing angle sin £ coming from neutrino masses, 
and it means that the /i-terms and VEVs have to be aligned, ui d //if oc v^jvj} '. In the MSSM case, a 
neutrino mass bound requires a strong alignment of the mixing term W ~ \x sin ^H U L^, 



sin£ < 3 x 10" Vl + tan 2 /3. 



;22) 



Therefore we simply neglect the mixing terms in the following analysis. 

We can rewrite the superpotential Eqs. (fTU|) ~ (fl3"l) in the completely same form with the physical 
fields, with new coupling constants defined as 



yd 

a 



A 



yd yd 
1 a J r b 



ye ye 

I b ) 1 b> 



V d yd 
U D 



\l yd 

d H> 

J D 



yd 
1 b' 



J D 



Vl \> y° 

U D 



Y° 



}^_ye , V L \ ye _ _ye 
d I b' ' A A > J c — I c i 



(23) 



ye 

1 b' 



J D 



7J A > A i - ~~d Y v 



Vd 



>D 



'D 



VL y d _ 

v d h 

D 



V d 



J D 



In what follows, we will discuss the phenomenology of these new superpotential written in unhatted 
fields. 



D. Fermion mass matrices and diagonalization 

We assume that VEVs take the form Eq. (fl4l) and (fl9l) . from which we obtain the fermion mass 
matrices. 

1. Quark sector 



The quark mass matrices are given by 



m 



m 



/ 

V2Y c u v%e- ld 3 
\ -Y^v u D e~ ieu 

( 

V2Y c d v$e- ie i 
-Y$v d D e- ied 



V2Y c u v%e~ i6 l 


Y$v u D e- ieu 

V2Y c d - d 




I b u D t 

Y h u v%e 



V2Y n u v%e- ie 3 



Y£v d D e 
Y^v d D e 



48 a 
-i0 d 



V2Y a » 



V 3 e 



(24) 



(25) 



In Ref 



22j . possibilities of the alignment have been discussed in the framework of U(l) horizontal symmetries. 
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We can bring the mass matrices above to the form 



/ q u /y u \ 



m 



m t 



\ 



-qu/yu o b u 
o K vlj 



(26) 



by 7r/4 rotations on the Qq doublet fermions and appropriate phase rotations. This mass matrix can 
be then diagonalized by orthogonal matrices 0\ R as 4 



' m u ' 

m c 
\ m t J 



(27) 



and similarly for m d . 

For the set of the parameters 



0d_ d d_ e u + u = _i.25,g„ 



0.0002150,6,, = 0.04440,6' = 0.09300, 



Vu 



0.99741, q d = 0.005040, b d = 0.02500, b' d = 0.7781, y d = 0.7970, 



:28i 



we obtain 

m u /m t 

| VcKM | 



1.11 x 10" 5 , m c /m t = 4.14 x 10" 
/ 0.97440 0.2267 0.00388 \ 

0.2265 0.9731 0.0421 
^ 0.00924 0.0412 0.9991 j 



,m d /m b = 1.22 x 10 -3 , m s /m b = 2.07 x 10" 



sin 2/9(^1 



0.723. 



(29) 



These values are consistent with present experimental values [8J. So, we see that the model can 
well reproduce the experimentally measured parameters. Moreover, since the CKM parameters 
and the quark masses are related to each other because of the family symmetry, we find that nine 
independent parameters (Eq. (f2"8"]l ) of the model can well describe ten physical observables: there is 
one prediction. An example of the prediction is |V^/V^|, whose experimental value has been obtained 
from the measurement of the mass difference Aitir of the B° meson [9[: 



Model : \V td /V ts \ 
Exp. : \V td /V ts \ 



0.21-0.23, 

+0.001 . . +0.008 . , 
0.208 exp. theo. 

-0.002 v ; -0.006 v 



(30) 



The form of the mass matrix is known as the next-neighbor interaction form 
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Finally, we give the unitary matrices that rotate the quarks for the choice of the parameters given 
in Eq. 



/ 0.706 0.0366 1.42 x 10~ 5 \ 



uL 



J 



+ e 



2iA8 u 



0.706 -0.0366 -1.42 x 10~ 5 


/ 0.0366 -0.705 0.0313 \ 

0.0366 -0.705 0.0313 

v -0.00231 e~ iA(>u 0.0441 e~ iA(>u 0.999 e~ iAeu j 



(31) 



U, 



dL 



( 0.695 0.130 0.00345 \ 
-0.695 -0.130 -0.00345 




+ e 



2iA8" 



0.130 
0.130 



/ 

-0.695 
-0.695 



0.0111 \ 
0.0111 



\ -0.00769 e- zAed 0.0146 e~ iAed 1.00 e~ iAei ) 



(32) 



U U R 



e 3 



0.0366 0.703 0.0657 \ 

0.0366 0.703 0.0657 

j 

( -0.706 



+ e 



\ 



0.706 
-0.00484 e~ iA0U 



0.0367 
-0.0367 



-6.73 x 10~ 6 \ 
6.74 x 10" 6 



-0.0928 e~ iA0U 0.996 e~ iA6u 



J 



(33) 



U, 



dR 



e 3 



_|_ e 2iA9 d +iQ$ 



( 0.134 0.427 0.548 \ 

0.134 0.427 0.548 
\ j 

( -0.675 
0.675 



0.212 -7.01xl0~ 5 \ 
-0.212 7.01 x 10~ 5 



(34) 



\ -0.232 e~ iAed -0.739 e~ iA(>d 0.633 e~ iAed J 
The unitary matrices above will be used when discussing R-parity violating processes in sect. III. 
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2. Lepton sector 



The mass matrix in the charged lepton sector is: 



m 



v D e 



V Y§ Y§ j 
It is diagonalized by the biunitary transformation: 



Ul L m e U eR 



m e 


0m 



\ 







One finds that U e z and U e R can be approximately written as 



17, 



ei 



-e P (l + e 2 
1 -e 



eV.- ■ 

2 
e 



:vV2)(i 



-e fj,J 

e 2 e 2 ) 



-y/2e 



1/V2 \ 
1/V2 
V2e e el J 



U, 



eR 



4/2) 



4/2 



V 



^(1-4: 

4^ 



4/2 y 



and small parameters e e , e„ are defined as 



3.42 x 10~ 3 , e 



— £ = 5.94 x 10 



-2 



In the limit m e = 0, the unitary matrix U e L becomes 

/ l/y/2 l/y/2~\ 

-l/y/2 l/y/2 

1 



(35) 



(36) 



(37) 



(3* 



(39) 



which is the origin for a maximal mixing of the atmospheric neutrinos. 

As for the neutrino sector, we assume that a see-saw mechanism [h]] takes place. However, we 
do not present the details of the neutrino sector here because there is no need to know it in the 
following analysis. We obtain some specific predictions of our model: (i) only an inverted mass 



hierarchy m U3 < m Ul , n 
(e, 3) element of the MNS matrix is given by \U e ^ 



is consistent with the experimental constraint lAm^l 



See Refs. 



< l A ^2 3 | ; 



%) the 



for details. 
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E. Soft supersymmetry breaking sector 



Next we consider the soft SUSY breaking sector, which respects the family symmetry. If three 
generations of a family is put into a one- dimensional and two-dimensional irreps of any dihedral 
group, then the soft scalar mass matrix for sfermions has always a diagonal form: 



/ 



m 



(q,£)LL 



rn 



\ 



(q,£)lL 










m 



(q,£)lL 










\ 



m 



(q,l)3L J 





/ 


m llR 








\ 









m llR 









V 








m lsR 


J 



(a = u,d,e). (40) 



Note that the mass of the first two generations are degenerated because of the family symmetry. 5 
Further, since the trilinear interactions (A-terms) are also Q§ invariant, the left-right mass matrices 
have the form 



( m Lfl)i,- 



u, d, e), 



(41) 



where Afs are free parameters of dimension one, and the fermion masses m's are given in Eq.(J2] 
(j25p and fl35|) . They are real, because we impose CP invariance at the Lagrangian level. 

We approximate that the squark and slepton masses are given only by Eq. (140j) . that is, trilinear 
terms Eq.( l4TI) are negligible j^. 



III. R-PARITY VIOLATION 



Since Q§ family symmetry controls the whole flavor structure of the model, the form of the R- 
parity violating couplings are also constrained by the family symmetry. We find that only possible 
trilinear couplings allowed by the symmetry can be written as 



W fi = XLsLjES + \' 1 L I (ia 2 ) IJ Q 3 D c J + X'Ma^uQjD 



(42) 



in the physical lepton doublet Lj defined in Eq. fllSp with the coupling constants in Eq. fl23j) . Here 
the superpotential is written in the flavor eigenstates, so the mixing matrices Eqs. (l3T!) ^ (!34|) . 
( |37|) and (|38|) should appear when we rotate the fermion components into their mass eigenstates. On 
the other hand, these matrices do not appear from sfermion components, because we approximate 
that sfermions are in the mass eigenstate basis. In the present model, there are only three R-parity 
violating trilinear interactions allowed by the family symmetry, and baryon number violating terms 
X'-j k U^DjDl are forbidden by the symmetry. It should be compared with the MSSM case in which 



5 SUSY flavor and CP problem can be avoided thanks to such partially degenerated (Eq. (JTDJ)) and aligned (Eq. pTj) ) 
soft terms because of the family symmetry [l|, Q. 



11 



there are 45 trilinear couplings. These interactions can generate a lot of new processes which have not 
been observed yet such as lepton flavor violating (LFV) processes, or new contributions to already 
observed processes. Many authors have studied phenomenology of R-parity violation and obtained 
constraints on each coupling constant corresponding to each process in the MSSM case. 6 

In this section, we obtain constraints on the coupling constants A, X[ 2 at the weak scale. Since the 
various new processes generated by the interactions Wr depend only on the three coupling constants, 
we can predict ratios of new processes independent of As. We will also find the ratios of the LFV 
processes in this section. As mentioned before, we assume that the R-parity violating couplings A 
and X[ 2 are real and positive 1 . 



A. Constraint on A 



In this subsection, we consider the constraint on XL 3 LjEj operator. The most stringent constraint 
on A is obtained from both fi — > eee and neutrinoless double beta decay, both the processes give 
similar bound. 

First we show the bound on A from /i — > eee process 0, HE 16, 17, 35]. The decay process 
£ m — > i^iji^ is generated by tree-level t- and u- channel sneutrino exchange (FIG. [[]), and its 
effective Lagrangian is given by 



C eff = \ 2 A L (liPdm) (IjPdh) + A 2 A R (liPd m ) ft^4) + (i <- j), 

where the coefficients are 





1 


A L = 














1 


A R = 













m 



£3L 



m hL 



IKk, 



IKm, 



(43) 

(44) 
(45) 



with mixing matrices in Eq. (1371) . (|38j) . In our approximation, sneutrino mass is the same as that of 
the left-handed slepton defined in Eq. fjlTJ]) . From this Lagrangian, the branching ratio of fi — > eee is 
given by, 

4A 4 



BR([i — > eee) 



6AG 2 F 



A L \ 2 + \A R \ 2 ] BR(n^ev e v^) 



(46) 



The requirement that this branching ratio should not exceed the experimental bound BR(fi 
eee) exp < 1.0 x 1CT 12 provides a constraint on A 



A < 1.4 x 10 



- 2 ( 



m 



VlOOGeV/ 



(47) 



See Ref. 
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and references therein. 

7 CP violation induced by the R-parity violating trilinear couplings in the soft SUSY breaking sector has been studied 
in Ref. 



39] 
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where we have assumed m 



eiL 



rrig 3L = mi L in order to forbid the contribution to FCNC processes 



from the soft scalar mass terms. 

Next we show the constraint from neutrinoless double beta d ecay . The coupling A can induce 
radiative neutrino mass at the one loop level shown in FIG. [T] [14, I20 ]: 



A 2 

(Mjy = -2j-^m ek [{U eL ) 3j (U eL ) Jk - (U eL ) 3 k(UeL)jj] x 



X 



(Ul R )ki(U eL h i (mj LR ) IJ F(m 2 ek , mj 1R , mj 1L ) - (Ul R ) kK (U eL ) Ki (m 2 SLR ) 3J F(m 2 ek , mj 1R , m\ L , 



where the loop function F is 

xlnx ylriy zlnz 

F{x,y,z) = - + 



(x-y)(x-z) (y-z)(y-x) (z - x)(z - y) 
\n(y/z) 



y-z 



as x -> 0. (49) 



We assume that the trilinear terms in the soft SUSY breaking sector are completely aligned with the 
Yukawa matrices, that is, 

mftfl ^ m e i e , (50) 

in which A e stands for the sum of the A-terms in Eq. (141 j) and //-terms which are not explicitly shown, 
and m e is the mass matrix of charged leptons Eq. (l35l) . The requirement that the effective neutrino 
mass (Ai u ) ee does not exceed a constraint from neutrinoless double beta decay provides a bound on 
A: 

A < 1.1 x 10- , M = A e ^° GeV l 2 m 4- (51) 

V0.35e^; \ M J rnl R -m% m 2 
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This is the most stringent constraint on A in the present model: other processes give weaker bounds. 
For example, a constraint from [i — - > is jl^] 

A < 0.22, (52) 

when slepton masses are lOOGeV. 



B. Constraints on \[ and X' 2 



Constraints on A' : 2 are obtained from neutral meson mixings 24|, |25|. 



on the products \\[ and AA 2 from leptonic decays of neutral mesons [26, 



26, 27, 28 



30 



31 



29, 35 , and 



32 



35 . Since 



both processes are generated at tree level, these give the most stringent bounds on \[ 2 . Although 
\i — e conversion in nuclei [l?], isj is also generated at tree level by \[ 2 , bounds from this process 
are weaker than those from neutral meson system. Therefore, we show the calculations only of the 
neutral meson system. 

The neutral meson mixing is generated at the tree level through the exchange of a sneutrino in 
both s- and t- channels (FIG. ED. For K° - K° mixing, the effective Hamiltonian is obtained as 



Ti 



A' A'* 



{d R s L )(d L s R ) 



(53) 



where 



A', 



ijk 



K{Ul R )ij{UdL)z k {icj 2 )ij + K(U ] dR ) r AU d L) 



IJ- 



(54) 



We require that these additional contributions to the mass difference of neutral K meson are smaller 
than its experimental value: 



Am 



n 

K° 



\A' I21 A 



I12\ 



S K om K ofKB A (m K ») < Am^o 



2 ™L 



where m K o and fx denote mass and the decay constant of the neutral K meson, and 

2 

/ m. rsn \ 

c 



m s {n) +m d (n) 



^ , B 4 (m K o) = 1.03, Am^T = (0.5292 ± 0.0009) x lO^ps -1 



with n = 2 GeV 28]. Then [3 



|A/2iA? 12 | < 4.5 x 10 



m 



I'L 



VIOOGeV/ 



The assumptions 8f — Of = and X[ = A' 2 lead to the most stringent constraints on X[ 2 : 



A' 1 = A 2 <3.1xlO- 3 ( m jL L r ) 
1 2 VIOOGeV/ 



(55) 



(56) 



(57) 



(5f 
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A' p A'(A) 




d% {in 



dj. 



X' 



di (£ rj 



dj(if) di- 



X' 



dj(£ f ) 



FIG. 2: R-parity violating contributions to the neutral meson mixings (or, leptonic dcays of the neutral 
mesons) in the s- and t- channels. The quantities in parenthesis represent the case of leptonic decays. 



Next, we consider the leptonic decays of the neutral K mesons: Kl,s - > e"e + ,/z~/i + or e T fi ± . At 
the quark level, these processes are interpreted as a transformation of a down-type quark-antiquark 
pair (di and dj) into a charged lepton-antilepton pair (if and i n ), which are shown in FIG. [2j The 
processes in the s- and t- channels are mediated by sneutrino and w-squark, respectively. In general, 



the effective Lagrangian for the process didj 

1 



C 



eff 



tfl n is written as 

Ajfn(djRl^d iR ) (IfLl/nL) 



32] 



+ 



~ — 2 ~ [(.djLdiR)(£fR£ n L)Bijfn + (djRdiL)(£fLi n R)Bji n f] 

I m M 



(59) 



where each coefficient is defined as 



A 



ijfn 



m 



+ 



A 



m 



- {(U eL ) Ln (ia 2 ) LK (U\ R ) 

■(UdR)3i(Ul R )j3(Ul L )f I (U eL ) In , 



(60) 



qlL 



Oijfn 



m 



M 



A 



m di + m dj m 



2 (ui R ) fK (u eL ) 3n AK^^ll)^, 



KJ 



eiL 



(61) 



and vtiM stands for the mass of the meson composed of quark(c/j)-antiquark(<ij) pair with mass 
and rridj- From the above effective Lagrangian, the decay rate of didj — > ifi n is given by 



T(didj 



"f n ) 



P I 

' !/ —CijfnJ {mli ~ ( m f ~ m n) 2 } {m 2 M - (m f + m n ) 2 }, 



2567rm 



M 



(62) 
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where 



Cijfn = (m 2 M - mj) \A ijfn m f + B ijfn \ 2 + (m 2 M - m 2 ) \A ijfn m n + B* inf \ 2 

| &ijfn + ^jinf | — | •Aijf n m n — Bijf n | 
|2 



| B ijfn m n - B* jinf m f \ + m f m n 



- \A ijfn m f -B* inf \ + \(m f + m n )Ai jf n\ J , (63) 

and J'm is the decay constant of the meson under consideration and m/ jn are the corresponding 
charged lepton masses. In the case of Kl decay, Aijf n should be replaced to (Aijf n — A/i/n) / V% with 
i — l,j — 2 , and similar for B{jf n and B* in j. 

We require that branching ratios should not exceed the experimental bound BR(Kl — > yU =F e ± ) exp < 
4.7 x 10 -12 and the experimantal value BR(Kl — > e~e + ) exp = 9 x 10~ 12 , therefore we obtain 

M' < 5.4xl0- 7 f m '" L V (64) 

1 VIOOGeV/ V ; 

from K L — ► /U^e 1 * 1 , and 

AA * < 11 x 10 " (k) 2 < 65 > 

from — > e~e + . The contributions including up type squarks, that is, Aijf n , are negligible in 
Kl — > / u T e ± because these terms only provide weaker constraints, and vanish in Kl — > e~e + because 
of the replacement mentioned above. 

For comparison, bounds on AA' X 2 from // — e conversion in nuclei are AA' : < 10~ 6 and AA' 2 < 10~ 7 , 
which are weaker than Eqs. (j64j) and (j65p . 



C. Predictions for Lepton Flavor Violating processes 

Although many processes can be generated by the R-parity violating interactions, we focus on the 
LFV decays £~ — ► tjtji^, where m = \i or r, in this subsection. As mentioned in the subsection 
A., the operator \L 3 L I E'f generates the decays £~ — * £^£j££ at tree level when A / 0. The 
other two operators in Eq. fj42|) . X' 12 LQD C , also generate the similar decay processes at one loop 
level through photon penguin diagrams shown in FIG. [31 but we found that the bounds on X[ 2 are 
stronger than that on A in the previous subsections. So we neglect contributions from X[ 2 operators 
to the decays £ m — > i^iji^- Moreover, flavor changing Z boson decay Z — > t~t^ induced by the 
R-parity violating bilinear terms can contribute to £~ — > £^£j£^ processes. Since branching ratios 



of these decays are propotional to sin 2 £, their effects are also negligible |23[]. Besides these R-parity 
violating contributions, there are two other contributions to these processes by Higgs bosons. Since 
the charged leptons couple to the neutral Higgs bosons, these Yukawa interactions generate LFV 



processes at one loop level 40] . However, these effects are enhanced only when tan/3 is large. So, 



16 



we assume that these are negligible because tan/3 is small enough. Moreover, since there are three 
generations of both up and down type Higgs doublet in this model, LFV processes mediated by the 
neutral Higgs bosons are generated at tree level. However, the branching ratio of the fi — > eee from 
these effects is BR ~ 10~ 16 because of the smallness of the Yukawa couplings when the neutral Higgs 
boson mass is lOOGeV. So, these contributions can also be negligible compared to those from XLLE C 
couplings unless A < 10~ 3 . Therefore, we can approximate that £ m — » £~£j£^ processes are induced 
at tree level only by A. In this approximation, the ratios of these processes are independent of A, but 
depend on the mixing matrices U e um which reflect the flavor structure of the model. Therefore we 



find some predictions of LFV decays 



3 k 



in our model. 



From the branching ratio Eq. (146p . we can easily find the ratios of processes in the approximation 
that all scalar masses are equal: 8 



BR(t 



eee) 



BR(t 
BR(t 



BR(t 



fififi) 
fifie) 



BR(T — > fififi) 

BR(fi — > eee) 



eee) 



4e 2 



0.014, 



~7< 

,2 



1-^ + 23 



e 2 - 2e 2 



T, 



A* 5 



0.99, 



0.0093, 



(66) 
(67) 
(68) 



where small parameters e e>At are given in Eq. (l3"9~|) and 7) t (r r ) stand for the lifetime of the ft(r) lepton. 
Also, BR(t — > fifie) means BR(r~ — > /x~/i~e + ), and similar for the other processes. One can obtain 
the ratios of other combinations from the branching ratios listed below: 



BR(n - 


-> eee) 


BR(t - 


-> eee) 


BR(t - 


> fififi) 


BR(t - 


-> ee/i) 


BR(r - 


■> fifie) 


BR(t - 


-> fiee) 


BR(r - 


■> /xe/x) 



2 (! - 2^1 + ^1- 2(1 -e> 



2 -2 



oc r r 



.2 
'7' 



m 7il + 2 e H 3 l 



4 (1 + 4 - 2^ - 4^eJ) m £3L 
2e 2 ,m~ 4 < r 1 - 2 



o 2 2 -2 -2 

2e e e M m |lI( m |3i 



-4 r> 2 -2 -2 

s-* — 2e,,m~ m-. 



1 



2 2 -4 



-e 4 

4 m 



m 



2e 2 t + 6e 2 e 2 ! + 



£3L' 

2 * 



m- 



£3L 



44 + 



3 -S 

4 ^ 



—4 1 o 2 2 — ^ — z 

m- 3L + 2e e e M m- 1L m- 3L 



(69) 
(70) 
(71) 
(72) 
(73) 
(74) 
(75) 



8 From the conditions to suppress \i — > e + 7 process from the scalar mass terms, slepton masses are required to be 
degenerated with mass differences of order 1O _1 0, 0|. 
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FIG. 3: The photon penguin contributions to the decays l m — * l i £j ^=j- These contributions can be 
negligible compared to the tree level processes induced by the coupling A. 



where the common factor is not shown explicitly. 



IV. CONCLUSION 

We have considered the properties of R-parity violating operators in a SUSY model with non- 
Abelian discrete Qq family symmetry. The family symmetry can reduce the number of parameters 
in the Yukawa sector, and explain the fermion masses and mixings between generations. It can 
also reduce the number of R-parity violating couplings and determine the form of those. Only 
three trilinear couplings are allowed, and the baryon number violating operators are forbidden by 
the symmetry in our model. We derived upper bounds on these couplings: A < O(10~ 2 ), A^ 2 < 
O(10~ 4 ), and obtained the predictions on the ratios of the LFV decays £~ — > £^£j£^ which do not 
depend unknown parameters. The results reflect the properties of the family symmetry because these 
predictions contain the mixing matrices of the charged lepton sector which is written by masses of the 
charged leptons. Our predictions can be testable at future experiments because the superB factory 



4l| or LHC will have the sensitivity BR ~ 10~ (8 ^ 9) for LFV r decays. 
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